The genetic regulatory architecture and epigenomic basis for age-related changes in rattlesnake venom.

Developmental phenotypic changes can evolve under selection imposed by age- and size-related ecological differences. Many of these changes occur through programmed alterations to gene expression patterns, but the molecular mechanisms and gene-regulatory networks underlying these adaptive changes remain poorly understood. Many venomous snakes, including the eastern diamondback rattlesnake (Crotalus adamanteus), undergo correlated changes in diet and venom expression as snakes grow larger with age, providing models for identifying mechanisms of timed expression changes that underlie adaptive life history traits. By combining a highly contiguous, chromosome-level genome assembly with measures of expression, chromatin accessibility, and histone modifications, we identified cis-regulatory elements and trans-regulatory factors controlling venom ontogeny in the venom glands of C. adamanteus. Ontogenetic expression changes were significantly correlated with epigenomic changes within genes, immediately adjacent to genes (e.g., promoters), and more distant from genes (e.g., enhancers). We identified 37 candidate transcription factors (TFs), with the vast majority being up-regulated in adults. The ontogenetic change is largely driven by an increase in the expression of TFs associated with growth signaling, transcriptional activation, and circadian rhythm/biological timing systems in adults with corresponding epigenomic changes near the differentially expressed venom genes. However, both expression activation and repression contributed to the composition of both adult and juvenile venoms, demonstrating the complexity and potential evolvability of gene regulation for this trait. Overall, given that age-based trait variation is common across the tree of life, we provide a framework for understanding gene-regulatory-network-driven life-history evolution more broadly.

The macroevolutionary singularity of snakes.

Snakes and lizards (Squamata) represent a third of terrestrial vertebrates and exhibit spectacular innovations in locomotion, feeding, and sensory processing. However, the evolutionary drivers of this radiation remain poorly known. We infer potential causes and ultimate consequences of squamate macroevolution by combining individual-based natural history observations (>60,000 animals) with a comprehensive time-calibrated phylogeny that we anchored with genomic data (5400 loci) from 1018 species. Due to shifts in the dynamics of speciation and phenotypic evolution, snakes have transformed the trophic structure of animal communities through the recurrent origin and diversification of specialized predatory strategies. Squamate biodiversity reflects a legacy of singular events that occurred during the early history of snakes and reveals the impact of historical contingency on vertebrate biodiversity.

Commensal protists in reptiles display flexible host range and adaptation to ectothermic hosts.

IMPORTANCE: Environmental factors like climate change and captive breeding can impact the gut microbiota and host health. Therefore, conservation efforts for threatened species may benefit from understanding how these factors influence animal microbiomes. Parabasalid protists are members of the mammalian microbiota that can modulate the immune system and impact susceptibility to infections. However, little is known about parabasalids in reptiles. Here, we profile reptile-associated parabasalids in wild and captive reptiles and find that captivity has minimal impact on parabasalid prevalence or diversity. However, because reptiles are cold-blooded (ectothermic), their microbiotas experience wider temperature fluctuation than microbes in warm-blooded animals. To investigate whether extreme weather patterns affect parabasalid-host interactions, we analyzed the gene expression in reptile-associated parabasalids and found that temperature differences significantly alter genes associated with host health. These results expand our understanding of parabasalids in this vulnerable vertebrate group and highlight important factors to be taken into consideration for conservation efforts.

Commensal protists in reptiles display flexible host range and adaptation to ectothermic hosts.

Parabasalid protists recently emerged as keystone members of the mammalian microbiota with important effects on their host's health. However, the prevalence and diversity of parabasalids in wild reptiles and the consequences of captivity and other environmental factors on these symbiotic protists are unknown. Reptiles are ectothermic, and their microbiomes are subject to temperature fluctuations, such as those driven by climate change. Thus, conservation efforts for threatened reptile species may benefit from understanding how shifts in temperature and captive breeding influence the microbiota, including parabasalids, to impact host fitness and disease susceptibility. Here, we surveyed intestinal parabasalids in a cohort of wild reptiles across three continents and compared these to captive animals. Reptiles harbor surprisingly few species of parabasalids compared to mammals, but these protists exhibited a flexible host-range, suggesting specific adaptations to reptilian social structures and microbiota transmission. Furthermore, reptile-associated parabasalids are adapted to wide temperature ranges, although colder temperatures significantly altered the protist transcriptomes, with increased expression of genes associated with detrimental interactions with the host. Our findings establish that parabasalids are widely distributed in the microbiota of wild and captive reptiles and highlight how these protists respond to temperature swings encountered in their ectothermic hosts.

Gut microbial ecology of Philippine gekkonids: ecoevolutionary effects on microbiome compositions.

Given the rapidly changing landscapes of habitats across the globe, a sound understanding of host-associated microbial communities and the ecoevolutionary forces that shape them is needed to assess general organismal adaptability. Knowledge of the symbiotic endogenous microbiomes of most reptilian species worldwide remains limited. We sampled gut microbiomes of geckos spanning nine species and four genera in the Philippines to (i) provide baseline data on gut microbiota in these host species, (ii) test for significant associations between host phylogenetic relationships and observed microbial assemblages, potentially indicative of phylosymbiosis, and (iii) identify correlations between multiple ecoevolutionary factors (e.g. species identity, habitat tendencies, range extents, and maximum body sizes) and gut microbiomes in Philippine gekkonids. We recovered no significant association between interspecific host genetic distances and observed gut microbiomes, providing limited evidence for phylosymbiosis in this group. Philippine gekkonid microbiomes were associated most heavily with host species identity, though marked variation among conspecifics at distinct sampling sites indicates that host locality influences gut microbiomes as well. Interestingly, individuals grouped as widespread and microendemic regardless of host species identity displayed significant differences in alpha and beta diversity metrics examined, likely driven by differences in rare OTU presence between groups. These results provide much needed insight in host-associated microbiomes in wild reptiles and the ecoevolutionary forces that structure such communities.

Distribution of snake species and snakebites in hotspots of Ethiopia.

INTRODUCTION: In Sub-Saharan Africa, snakebites are a public health problem. In Ethiopia, clinical cases have been described, but little information exists on snakebites burden and its geographical distribution. The aim of this study was to document the spatial distribution of venomous snakes and snakebites in Ethiopia. METHODOLOGY: In a cross-sectional observational study, venomous snakes were collected during snake catching activities in six Ethiopian hotspot areas between April 2015 and September 2020. Species and habitat were described. In the hotspot areas, routine health information data on reported snakebites was collected in 78 districts and subsequently used to map annual incidence per district. RESULTS: A total of 333 snakes were collected and 14 species were identified. The most prevalent species were Bitis arietans, Bitis arietanus somalica, Echis pyramidum, known as vipers, and Naja pallida, known as cobra. The highest number of snakes (75) was observed in the Northwest and Eastern parts of Ethiopia, mostly in cultivation and man-made farm land, wooded and moist dry savanna. In each hotspot a wide variety of species was observed, although composition was different. The highest snakebite incidence overlapped with the high snakes densities in Northwest Ethiopia. The snakebite annual average incidence at district level was very heterogeneous and ranged from < 15 cases/100,000 inhabitants (44% of the districts) to 309.2 cases/100,000 inhabitants. CONCLUSIONS: Snake diversity and distribution, linked to high incidence of snakebites in the hotspots, suggests a close interconnection between human, animal and environmental systems and could inform the need for antivenoms per geographical locality.

A comprehensive phylogeny of dwarf geckos of the genus Lygodactylus, with insights into their systematics and morphological variation.

The 71 currently known species of dwarf geckos of the genus Lygodactylus are a clade of biogeographic interest due to their occurrence in continental Africa, Madagascar, and South America. Furthermore, because many species are morphologically cryptic, our knowledge of species-level diversity within this genus is incomplete, as indicated by numerous unnamed genetic lineages revealed in previous molecular studies. Here we provide an extensive multigene phylogeny covering 56 of the named Lygodactylus species, four named subspecies, and 34 candidate species of which 19 are newly identified in this study. Phylogenetic analyses, based on ∼10.1 kbp concatenated sequences of eight nuclear-encoded and five mitochondrial gene fragments, confirm the monophyly of 14 Lygodactylus species groups, arranged in four major clades. We recover two clades splitting from basal nodes, one comprising exclusively Malagasy species groups, and the other containing three clades. In the latter, there is a clade with only Madagascar species, which is followed by a clade containing three African and one South American species groups, and its sister clade containing six African and two Malagasy species groups. Relationships among species groups within these latter clades remain weakly supported. We reconstruct a Lygodactylus timetree based on a novel fossil-dated phylotranscriptomic tree of squamates, in which we included data from two newly sequenced Lygodactylus transcriptomes. We estimate the crown diversification of Lygodactylus started at 46 mya, and the dispersal of Lygodactylus among the main landmasses in the Oligocene and Miocene, 35-22 mya, but emphasize the wide confidence intervals of these estimates. The phylogeny suggests an initial out-of-Madagascar dispersal as most parsimonious, but accounting for poorly resolved nodes, an out-of-Africa scenario may only require one extra dispersal step. More accurate inferences into the biogeographic history of these geckos will likely require broader sampling of related genera and phylogenomic approaches to provide better topological support. A survey of morphological characters revealed that most of the major clades and species groups within Lygodactylus cannot be unambiguously characterized by external morphology alone, neither by unique character states nor by a diagnostic combination of character states. Thus, any future taxonomic work will likely benefit from integrative, phylogenomic approaches.

Gut microbial ecology of the Critically Endangered Fijian crested iguana (Brachylophus vitiensis): Effects of captivity status and host reintroduction on endogenous microbiomes.

Animals often exhibit distinct microbial communities when maintained in captivity as compared to when in the wild. Such differentiation may be significant in headstart and reintroduction programs where individuals spend some time in captivity before release into native habitats. Using 16S rRNA gene sequencing, we (i) assessed differences in gut microbial communities between captive and wild Fijian crested iguanas (Brachylophus vitiensis) and (ii) resampled gut microbiota in captive iguanas released onto a native island to monitor microbiome restructuring in the wild. We used both cloacal swabs and fecal samples to further increase our understanding of gut microbial ecology in this IUCN Critically Endangered species. We found significant differentiation in gut microbial community composition and structure between captive and wild iguanas in both sampling schemes. Approximately two months postrelease, microbial communities in cloacal samples from formerly captive iguanas closely resembled wild counterparts. Interestingly, microbial communities in fecal samples from these individuals remained significantly distinct from wild conspecifics. Our results indicate that captive upbringings can lead to differences in microbial assemblages in headstart iguanas as compared to wild individuals even after host reintroduction into native conditions. This investigation highlights the necessity of continuous monitoring of reintroduced animals in the wild to ensure successful acclimatization and release.

Venomous Snakes Reveal Ecological and Phylogenetic Factors Influencing Variation in Gut and Oral Microbiomes.

The gastrointestinal tract (GIT) of vertebrates contains a series of organs beginning with the mouth and ending with the anus or cloacal opening. Each organ represents a unique environment for resident microorganisms. Due to their simple digestive anatomy, snakes are good models for studying microbiome variation along the GIT. Cloacal sampling captures the majority of the microbial diversity found in the GIT of snakes-yet little is known about the oral microbiota of snakes. Most research on the snake mouth and gut microbiota are limited to studies of a single species or captive-bred individuals. It therefore remains unclear how a host's life history, diet, or evolutionary history correlate with differences in the microbial composition within the mouths and guts of wild snakes. We sampled the mouth and gut microbial communities from three species of Asian venomous snakes and utilized 16S rRNA microbial inventories to test if host phylogenetic and ecological differences correlate with distinct microbial compositions within the two body sites. These species occupy three disparate habitat types: marine, semi-arboreal, and arboreal, our results suggest that the diversity of snake mouth and gut microbial communities correlate with differences in both host ecology and phylogeny.

Phylogenetically diverse diets favor more complex venoms in North American pitvipers.

The role of natural selection in the evolution of trait complexity can be characterized by testing hypothesized links between complex forms and their functions across species. Predatory venoms are composed of multiple proteins that collectively function to incapacitate prey. Venom complexity fluctuates over evolutionary timescales, with apparent increases and decreases in complexity, and yet the causes of this variation are unclear. We tested alternative hypotheses linking venom complexity and ecological sources of selection from diet in the largest clade of front-fanged venomous snakes in North America: the rattlesnakes, copperheads, cantils, and cottonmouths. We generated independent transcriptomic and proteomic measures of venom complexity and collated several natural history studies to quantify dietary variation. We then constructed genome-scale phylogenies for these snakes for comparative analyses. Strikingly, prey phylogenetic diversity was more strongly correlated to venom complexity than was overall prey species diversity, specifically implicating prey species' divergence, rather than the number of lineages alone, in the evolution of complexity. Prey phylogenetic diversity further predicted transcriptomic complexity of three of the four largest gene families in viper venom, showing that complexity evolution is a concerted response among many independent gene families. We suggest that the phylogenetic diversity of prey measures functionally relevant divergence in the targets of venom, a claim supported by sequence diversity in the coagulation cascade targets of venom. Our results support the general concept that the diversity of species in an ecological community is more important than their overall number in determining evolutionary patterns in predator trait complexity.

The Tiger Rattlesnake genome reveals a complex genotype underlying a simple venom phenotype.

Variation in gene regulation is ubiquitous, yet identifying the mechanisms producing such variation, especially for complex traits, is challenging. Snake venoms provide a model system for studying the phenotypic impacts of regulatory variation in complex traits because of their genetic tractability. Here, we sequence the genome of the Tiger Rattlesnake, which possesses the simplest and most toxic venom of any rattlesnake species, to determine whether the simple venom phenotype is the result of a simple genotype through gene loss or a complex genotype mediated through regulatory mechanisms. We generate the most contiguous snake-genome assembly to date and use this genome to show that gene loss, chromatin accessibility, and methylation levels all contribute to the production of the simplest, most toxic rattlesnake venom. We provide the most complete characterization of the venom gene-regulatory network to date and identify key mechanisms mediating phenotypic variation across a polygenic regulatory network.

Congruence and Conflict in the Higher-Level Phylogenetics of Squamate Reptiles: An Expanded Phylogenomic Perspective.

Genome-scale data have the potential to clarify phylogenetic relationships across the tree of life but have also revealed extensive gene tree conflict. This seeming paradox, whereby larger data sets both increase statistical confidence and uncover significant discordance, suggests that understanding sources of conflict is important for accurate reconstruction of evolutionary history. We explore this paradox in squamate reptiles, the vertebrate clade comprising lizards, snakes, and amphisbaenians. We collected an average of 5103 loci for 91 species of squamates that span higher-level diversity within the clade, which we augmented with publicly available sequences for an additional 17 taxa. Using a locus-by-locus approach, we evaluated support for alternative topologies at 17 contentious nodes in the phylogeny. We identified shared properties of conflicting loci, finding that rate and compositional heterogeneity drives discordance between gene trees and species tree and that conflicting loci rarely overlap across contentious nodes. Finally, by comparing our tests of nodal conflict to previous phylogenomic studies, we confidently resolve 9 of the 17 problematic nodes. We suggest this locus-by-locus and node-by-node approach can build consensus on which topological resolutions remain uncertain in phylogenomic studies of other contentious groups. [Anchored hybrid enrichment (AHE); gene tree conflict; molecular evolution; phylogenomic concordance; target capture; ultraconserved elements (UCE).].

Phylogenetically diverse diets favor more complex venoms in North American pitvipers

The role of natural selection in the evolution of trait complex-ity can be characterized by testing hypothesized links betweencomplex forms and their functions across species. Predatory ven-oms are composed of multiple proteins that collectively function toincapacitate prey. Venom complexity fluctuates over evolutionarytimescales, with apparent increases and decreases in complexity,and yet the causes of this variation are unclear. We tested alterna-tive hypotheses linking venom complexity and ecological sourcesof selection from diet in the largest clade of front-fanged ven-omous snakes in North America: the rattlesnakes, copperheads,cantils, and cottonmouths. We generated independent transcrip-tomic and proteomic measures of venom complexity and collatedseveral natural history studies to quantify dietary variation. Wethen constructed genome-scale phylogenies for these snakes forcomparative analyses. Strikingly, prey phylogenetic diversity wasmore strongly correlated to venom complexity than was overallprey species diversity, specifically implicating prey species’ diver-gence, rather than the number of lineages alone, in the evolutionof complexity. Prey phylogenetic diversity further predicted tran-scriptomic complexity of three of the four largest gene familiesin viper venom, showing that complexity evolution is a concertedresponse among many independent gene families. We suggest thatthe phylogenetic diversity of prey measures functionally relevantdivergence in the targets of venom, a claim supported by sequencediversity in the coagulation cascade targets of venom. Our resultssupport the general concept that the diversity of species in an eco-logical community is more important than their overall number indetermining evolutionary patterns in predator trait complexity.

Gradual and Discrete Ontogenetic Shifts in Rattlesnake Venom Composition and Assessment of Hormonal and Ecological Correlates.

Ontogenetic shifts in venom occur in many snakes but establishing their nature as gradual or discrete processes required additional study. We profiled shifts in venom expression from the neonate to adult sizes of two rattlesnake species, the eastern diamondback and the timber rattlesnake. We used serial sampling and venom chromatographic profiling to test if ontogenetic change occurs gradually or discretely. We found evidence for gradual shifts in overall venom composition in six of eight snakes, which sometimes spanned more than two years. Most chromatographic peaks shift gradually, but one quarter shift in a discrete fashion. Analysis of published diet data showed gradual shifts in overall diet composition across the range of body sizes attained by our eight study animals, while the shifts in abundance of different prey classes varied in form from gradual to discrete. Testosterone concentrations were correlated with the change in venom protein composition, but the relationship is not strong enough to suggest causation. Venom research employing simple juvenile versus adult size thresholds may be failing to account for continuous variation in venom composition lifespan. Our results imply that venom shifts represent adaptive matches to dietary shifts and highlight venom for studies of alternative gene regulatory mechanisms.

A new species of African Snake-Eyed Skink (Scincidae: Panaspis) from Ethiopia.

A recent molecular phylogenetic revision of the snake-eyed skinks (genus Panaspis Cope, 1868) uncovered extensive cryptic diversity, including several new species from throughout sub-Saharan Africa. Here, we describe one of these from Ethiopia as Panaspis annettesabinae sp. nov. This description is based on a previous molecular phylogeny and morphological, scalation, and coloration data collected from the type specimen. Phylogenetic analyses place the species alone in what we term the P. annettesabinae species group from Ethiopia. This group forms the sister lineage to a large southern African radiation and suggests a potential northern origin for much of the extant diversity of Panaspis. Many new taxa have recently been discovered in the genus and region, and there are several historical records of Panaspis from elsewhere in Ethiopia. Thus, we suggest that the range of this new species (known only from a single specimen at present) may be much larger, and that additional undescribed species may exist in northern sub-Saharan Africa.

Availability of Recently Described Ptychadena (Anura: Ptychadenidae) Nomina from Ethiopia.

Frogs of the genus Ptychadena (Boulenger, 1917) have long been identified as harboring cryptic diversity and comprising numerous species-complexes (Largen 1997; Zimkus et al. 2017), and many authors have recognized the particularly high hidden richness in the Ethiopian highlands (Largen 1997; see Largen Spawls 2010 and refs. within). This cryptic diversity was confirmed by recent molecular studies (Freilich et al. 2014; Smith et al. 2017a, Reyes-Velasco et al. 2018). Those authors identified a congruent set of evolutionarily distinct candidate species using both mitochondrial and nuclear DNA, and described the geographic and ecological isolation of these species in detail (Freilich et al. 2014; Smith et al. 2017a).

Phylogenetic and spatial distribution of evolutionary diversification, isolation, and threat in turtles and crocodilians (non-avian archosauromorphs).

BACKGROUND: The origin of turtles and crocodiles and their easily recognized body forms dates to the Triassic and Jurassic. Despite their long-term success, extant species diversity is low, and endangerment is extremely high compared to other terrestrial vertebrate groups, with ~ 65% of ~ 25 crocodilian and ~ 360 turtle species now threatened by exploitation and habitat loss. Here, we combine available molecular and morphological evidence with statistical and machine learning algorithms to present a phylogenetically informed, comprehensive assessment of diversification, threat status, and evolutionary distinctiveness of all extant species. RESULTS: In contrast to other terrestrial vertebrates and their own diversity in the fossil record, the recent extant lineages of turtles and crocodilians have not experienced any global mass extinctions or lineage-wide shifts in diversification rate or body-size evolution over time. We predict threat statuses for 114 as-yet unassessed or data-deficient species and identify a concentration of threatened turtles and crocodilians in South and Southeast Asia, western Africa, and the eastern Amazon. We find that unlike other terrestrial vertebrate groups, extinction risk increases with evolutionary distinctiveness: a disproportionate amount of phylogenetic diversity is concentrated in evolutionarily isolated, at-risk taxa, particularly those with small geographic ranges. Our findings highlight the important role of geographic determinants of extinction risk, particularly those resulting from anthropogenic habitat-disturbance, which affect species across body sizes and ecologies. CONCLUSIONS: Extant turtles and crocodilians maintain unique, conserved morphologies which make them globally recognizable. Many species are threatened due to exploitation and global change. We use taxonomically complete, dated molecular phylogenies and various approaches to produce a comprehensive assessment of threat status and evolutionary distinctiveness of both groups. Neither group exhibits significant overall shifts in diversification rate or body-size evolution, or any signature of global mass extinctions in recent, extant lineages. However, the most evolutionarily distinct species tend to be the most threatened, and species richness and extinction risk are centered in areas of high anthropogenic disturbance, particularly South and Southeast Asia. Range size is the strongest predictor of threat, and a disproportionate amount of evolutionary diversity is at risk of imminent extinction.

Sexual Dichromatism Drives Diversification within a Major Radiation of African Amphibians.

Theory predicts that sexually dimorphic traits under strong sexual selection, particularly those involved with intersexual signaling, can accelerate speciation and produce bursts of diversification. Sexual dichromatism (sexual dimorphism in color) is widely used as a proxy for sexual selection and is associated with rapid diversification in several animal groups, yet studies using phylogenetic comparative methods to explicitly test for an association between sexual dichromatism and diversification have produced conflicting results. Sexual dichromatism is rare in frogs, but it is both striking and prevalent in African reed frogs, a major component of the diverse frog radiation termed Afrobatrachia. In contrast to most other vertebrates, reed frogs display female-biased dichromatism in which females undergo color transformation, often resulting in more ornate coloration in females than in males. We produce a robust phylogeny of Afrobatrachia to investigate the evolutionary origins of sexual dichromatism in this radiation and examine whether the presence of dichromatism is associated with increased rates of net diversification. We find that sexual dichromatism evolved once within hyperoliids and was followed by numerous independent reversals to monochromatism. We detect significant diversification rate heterogeneity in Afrobatrachia and find that sexually dichromatic lineages have double the average net diversification rate of monochromatic lineages. By conducting trait simulations on our empirical phylogeny, we demonstrate that our inference of trait-dependent diversification is robust. Although sexual dichromatism in hyperoliid frogs is linked to their rapid diversification and supports macroevolutionary predictions of speciation by sexual selection, the function of dichromatism in reed frogs remains unclear. We propose that reed frogs are a compelling system for studying the roles of natural and sexual selection on the evolution of sexual dichromatism across micro- and macroevolutionary timescales.

A species-level phylogeny of Trachylepis (Scincidae: Mabuyinae) provides insight into their reproductive mode evolution.

Trachylepis (Mabuyinae) includes ∼80 species of fully-limbed skinks found primarily in Africa and Madagascar, but a robust species-level phylogeny for this genus is lacking and this impedes studies on a wide-range of topics from biogeography to character evolution. Trachylepis and its close relatives (which together form the Mabuya group or Mabuyinae) are notable in that they have undergone multiple transitions and remarkable specializations in their reproductive modes. A Trachylepis phylogeny will be particularly useful for investigating reproductive evolution, because it includes species that exhibit oviparity, viviparity, and bimodal parity (species with both oviparous and viviparous populations). We sequenced DNA at four mitochondrial and five nuclear loci for 67 (∼84% of) Trachylepis species to infer a phylogeny for this genus. We performed stochastic character mapping of parity mode under a variety of parity mode transition models to infer ancestral parity mode states and the number and type of parity mode transitions. We recovered a strongly supported phylogeny of Trachylepis that is generally consistent with earlier phylogenetic studies. The best-fit model of reproductive mode evolution supports an oviparous ancestor for Trachylepis, and supports at least three viviparity to oviparity transitions. We compared parity mode evolution under the overall best-fit model (no constraints on parity mode transitions) to the best-fit model among the subset of models that assume viviparity to oviparity transitions are impossible. Our results support a model of reproductive evolution that allows for reversibility from viviparity to oviparity, a process that is not generally accepted. Alternatively, the best-fit model of evolution among the set of models that eliminate reversals from viviparity to oviparity suggests that bimodal reproduction may have persisted for millions of years within multiple lineages.

The emerging field of venom-microbiomics for exploring venom as a microenvironment, and the corresponding Initiative for Venom Associated Microbes and Parasites (iVAMP).

Venom is a known source of novel antimicrobial natural products. The substantial, increasing number of these discoveries have unintentionally culminated in the misconception that venom and venom-producing glands are largely sterile environments. Culture-dependent and -independent studies on the microbial communities in venom microenvironments reveal the presence of archaea, algae, bacteria, fungi, protozoa, and viruses. Venom-centric microbiome studies are relatively sparse to date with the adaptive advantages that venom-associated microbes might offer to their hosts, or that hosts might provide to venom-associated microbes, remaining largely unknown. We highlight the potential for the discovery of venom microbiomes within the adaptive landscape of venom systems. The considerable number of convergently evolved venomous animals, juxtaposed with the comparatively few known studies to identify microbial communities in venom, provides new possibilities for both biodiversity and therapeutic discoveries. We present an evidence-based argument for integrating microbiology as part of venomics (i.e., venom-microbiomics) and introduce iVAMP, the Initiative for Venom Associated Microbes and Parasites (https://ivamp-consortium.github.io/), as a growing collaborative consortium. We express commitment to the diversity, inclusion and scientific collaboration among researchers interested in this emerging subdiscipline through expansion of the iVAMP consortium.